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O ! ABSTRACT 

^ I We present a multi-wavelength study of NGC 4330, a highly-inclined spiral galaxy in the Virgo Cluster which 

^ . is a clear example of strong, ongoing ICM-ISM ram pressure stripping. The HI has been removed from well 

■ within the undisturbed old stellar disk, to 50% - 65 % of R25 . Multi-wavelength data (WI ¥1^0 B VR-Ha, VLAQ 

21 -cm HI and radio continuum, and GALEXQ NUV and FUV) reveal several one-sided extraplanar features 
CS| ' likely caused by ram pressure at an intermediate disk-wind angle. At the leading edge of the interaction, the 

' Ha and dust extinction curve sharply out of the disk in a remarkable and distinctive "upturn" feature that may 

! be generally useful as a diagnostic indicator of active ram pressure. On the trailing side, the ISM is stretched 

i out in a long tail which contains 10% of the galaxy's total HI emission, 6 - 9% of its NUV-FUV emission, but 

' only 2% of the Ha. The centroid of the HI tail is downwind of the UV/Ha tail, suggesting that the ICM wind 

I has shifted most of the ISM downwind over the course of the past 10 - 300 Myr. Along the major axis, the 

disk is highly asymmetric in the UV, but more symmetric in Ha and HI, also implying recent changes in the 
distributions of gas and star formation. The UV-optical colors indicate very different star formation histories 
for the leading and trailing sides of the galaxy. On the leading side, a strong gradient in the UV-optical colors 
of the gas-stripped disk suggests that it has taken 200-400 Myr to strip the gas from a radius of >8 to 5 kpc, 
but on the trailing side there is no age gradient. All our data suggest a scenario in which NGC 4330 is falling 
into cluster center for first time and has experienced a significant increase in ram pressure over the last 200-400 
^ , Myr 

\—i ' Subject headings: galaxies: individual (NGC 4330) — galaxies: spiral — galaxies: evolution — galaxies: ISM 

5^ I — galaxies: structure — galaxies: interactions 



1. EVITRODUCTION 

There is ample evidence that cluster galaxies are strongly 
affected by their environment over time. Denser environ- 
ments in the nearby universe have a higher fraction of ellip- 
tical and SO galaxies (the density-morphology relationship of 
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Dressier 1980) and the spirals are systematically Hl-deficient 
(Haynes & Giovanelli 1986, Solanes et al. 2001), clusters 
at z>0.1 have a higher fraction of blue, star-forming galax- 
ies than nearby clusters (Butcher & Oemler 1978, 1984), and 
at z^O.5, the fraction of spirals increases while the fraction of 
SO's decreases relative to that in the local universe (Dressier et 
al. 1997). This suggests that one or more cluster processes act 
to transform spiral galaxies into earlier-type spirals and SO's. 
Many effects have been presented as possible contributors, in- 
cluding galaxy-galaxy harrassment (Bekki 1999, Moore et al. 
1998), galaxy-galaxy collisions (Kenney et al. 2008), galaxy 
interactions with the cluster potential (Moore et al. 1999), 
ram pressure stripping (Gunn & Gott 1972, Abadi et al. 1999, 
Poggianti et al. 1999), turbulent viscous stripping (Nulsen, 



2 



Abramson et al. 



1982), and starvation (Larson et al. 1980). There remains 
significant debate about the relative importance of these pro- 
cesses. 

Several recent studies of clusters out to a redshift of z^O.8 
support the importance of ICM-ISM interactions as a driver 
of cluster galaxy evolution. Many studies have found cluster 
populations of disky, quiescent galaxies with structure similar 
to spiral galaxies but with little or no star formation (e.g., van 
den Bergh 1976, Balogh et al. 1998, van der Wei et al. 2010). 
Such galaxies are consistent with gas stripping. Moran et al. 
(2007) present evidence from two z^O.5 clusters that there 
are multiple pathways for star formation to be quenched in 
spirals, depending on the properties of the cluster, with ICM- 
ISM stripping emerging as a major contributing factor in sit- 
uations where the ICM is dense enough for stripping to be 
effective. From a spectral analysis of galaxies in various en- 
vironments at z = 0.4 - 0.8, Poggianti et al. (2009) find that 
that clusters are the preferred location of k-na galaxies. Such 
galaxies have had their star formation was quenched between 
5x 10^ and 1.5 x 10^ years prior to observation, making them 
a likely transition phase between actively star-forming spirals 
and passively evolving SO and Sa galaxies. Since the star for- 
mation quenching efficiency increases with cluster velocity 
dispersion, Poggianti et al. (2009) propose this transforma- 
tion is likely the result of ICM-lSM interactions. These stud- 
ies document how key galaxy properties differ for large sam- 
ples of cluster galaxies, but do not by themselves offer direct 
evidence of the mechanism transforming cluster galaxies. 

In order to understand the impact of ram pressure stripping 
on galaxy evolution, we need to be able to identify galaxies 
which are clearly experiencing stripping or did at some point 
in the past, and determine which of a given galaxy's char- 
acteristics are due to ram pressure stripping. In particular, we 
are interested in finding out how the multi-phase ISM behaves 
during stripping, how long it takes to strip a galaxy, and how 
ram pressure and the associated stripping affects the rate and 
distribution of star formation. Through detailed studies of the 
best cases of ram pressure stripping in the nearby Virgo clus- 
ter, we hope to identify key parameters of the ram pressure 
stripping interaction, such as the current strength of ram pres- 
sure, the angle between the ICM wind and the disk, the time 
evolution of these quantities, and the evolutionary stage of the 
stripping event. If we can quantify key stripping parameters 
for enough galaxies, we can understand the impact of strip- 
ping on galaxy evolution in different environments. 

There have been many ram pressure stripping simulations 
over the past few years (e.g., Quilis et al. 2000, Schulz 
& Struck 2001, Roediger & Bruggen 2006, Tonnesen et al. 
2007, Tonnesen & Bryan 2008, 2009, 2010, Ki-onberger et 
al. 2008, Vollmer 2009, Jachym et al. 2009; also see re- 
view by E. Roediger 2009). They are all in general agree- 
ment that in most situations, ram pressure stripping is well- 
approximated by the simple Gunn & Gott (1972) formula, 
which fails only in cases of nearly edge-on stripping and short, 
impulsive episodes of ram pressure (Jachym et al. 2009). All 
of the simulations predict that the gas stripped from the disk 
will form a tail downwind of the galaxy, and all of the groups 
predict a lack of star formation in the stripped outer disks. 
Kronberger et al. (2008) and Tonnesen & Bryan (2009) both 
predict enhanced star formation rates due to ICM pressure in 
certain cases. Progress has been made in modeling the effects 
of the stripping process on the multi-phase ISM. For exam- 
ple, Tonnesen & Bryan (2009, 2010) model the effects of ram 



pressure on a clumpy ISM with self-gravity and cooling, and 
also account for the effects of heating. However, more work 
remains to be done. Many complicated physical processes 
relevant to the ISM, particularly star formation, are simulated 
using simple prescriptions, and observations are necessary to 
determine what actually happens. 

There have also been spectacular examples of individual 
galaxies being stripped in rich clusters between z^O.2 and 
0.5 (Owen et al. 2006, Cortese et al. 2007, Sun et al. 2007), 
but at such distances we lack the resolution to see what re- 
ally happens. The Virgo cluster is the nearest moderately rich 
cluster, and provides the best resolution to observe cluster pro- 
cesses at work. For this reason, we have carried out the VIVA 
(VLA Imaging of Virgo spirals in Atomic gas) survey, which 
presents high-quality HI imaging of 53 spirals in the Virgo 
cluster (Chung et al. 2009). Using data from the VIVA survey, 
Chung et al. (2007) noted the presence of HI tails in at least 
25% of large spiral galaxies, a total of 7 galaxies, near the 
cluster's virial radius, providing striking evidence of galaxies 
actively losing their gas as they plummet towards the cluster 
center. One of these galaxies was NGC 4330, and we have 
followed up with a detailed multi-wavelength study. 

In this paper, we show that NGC 4330 is one of the best 
examples of ram pressure stripping actively transforming a 
spiral galaxy in the nearby Virgo cluster. There is a wealth 
of evidence that it is being actively stripped, and several in- 
teresting phenomena can be clearly seen, some for the first 
time. In this multi-wavelength study, we examine the effects 
of stripping on the galaxy's current star formation, recent star 
formation history, and the stripping of the HI and dust compo- 
nents of its ISM. In Section 12 we introduce the galaxy NGC 
4330. In Section[3] we describe the observations taken in op- 
tical broad-band and Ha narrow-band filters, HI, and the UV. 
In SectionlH we describe the galaxy's stellar distribution and 
optical dust extinction in the BVR bands, the HI distribution 
and kinematics, and the Ha and UV distributions. In Section 
|5] we discuss the interpretation of the galaxy's properties in 
the context of an ongoing ICM-ISM interaction. In Section 
121 we summarize our findings. Throughout this paper, we as- 
sume a distance to the Virgo Cluster of 16 Mpc (Yasuda et al. 
1997). 

2. THE GALAXY NGC 4330 

NGC 4330 is a late-type, small-bulge, highly-inclined spi- 
ral with a maximum rotation speed of ~ 180 km/s, a lumi- 
nosity of ^ 0.3L*, and an HI mass of 4.5 x 10^ Mq (Chung 
et al. 2009). It is located at a projected distance of 2° (600 
kpc) from M87 (see Figure[T]), in a region near the Virgo clus- 
ter's virial radius where many galaxies are HI deficient and 
some have HI gas tails likely resulting from ICM ram pressure 
(Chung et al. 2007). NGC 4330 fulfills both of these criteria, 
making it a strong candidate for ongoing ram pressure strip- 
ping. This galaxy was initially noticed in GALEX UV images 
because of its UV-bright tail. Chung et al. (2007) find that this 
galaxy possesses an HI tail that extends to the south of the disk 
plane on the western side, pointing roughly away from M87. 
Using the assumptions of Gunn & Gott (1972), Chung et al. 
(2007) estimate that the ICM ram pressure on the galaxy at 
this location in the cluster is sufficient to cause the observed 
stripping of HI. The tail direction suggests that NGC 4330 
is currently falling in towards the cluster core, and the pres- 
ence of a dense gas tail strongly suggests active stripping. The 
galaxy's HI deficiency parameter is 0.8 (Chung et al. 2007), 
meaning that its HI mass is only 15% of that expected for a 
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field spiral of similar size. NGC 4330 has a heliocentric veloc- 
ity of 1569 km/s, which exceeds the Virgo Cluster's mean he- 
liocentric systemic velocity of ^11 00 km s^'. This indicates 
that NGC 4330 is likely moving roughly towards the core of 
the cluster from the front. With an inclination of 84°, NGC 
4330 is nearly edge-on. This makes it an excellent galaxy in 
which to examine the effects of ICM pressure on the galaxy's 
recent star formation and HI gas morphology, since we can 
distinguish between planar and extraplanar gas. 

NGC 4330 displays a local radio continuum deficit (Figure 
121, which is believed to be a tracer of ram pressure (Murphy 
et al. 2009). In a sample of Virgo spiral galaxies examined by 
Murphy et al. (2009), NGC 4330 had the sti'ongest local radio 
deficit, suggesting that it is experiencing strong, active ram 
pressure. Murphy et al. (2009) conclude that the radio deficit 
regions are approximately coincident with the galaxies' lead- 
ing edges, where the ICM is actively exerting force on the 
ISM. Since the galaxy's tail, which should be roughly aligned 
with the ICM wind direction, is on the side of the galaxy op- 
posite the radio deficit region, we feel confident calling the 
side of the galaxy with the radio deficit region "upwind" and 
the side with the tail "downwind" of the major axis. 

Gravitational interactions cannot account for the undis- 
turbed stellar disk (see Section 14. It and disturbed gas dis- 
tribution in NGC 4330. In any case, there are no extremely 
close neighbors that would be likely to cause strong gravita- 
tional interactions. The only Virgo galaxy within 25' (^100 
kpc projected) and 500 km s"' of NGC 4330 is VCC 0706, an 
early-type dwarf galaxy at a projected distance of 77 kpc and 
velocity difference of 75 km s^'. Since it is 4.2 magnitudes 
fainter than NGC 4330 (Br = 17.3, and NGC 4330 has Br = 
13.1), it is too small to significantly perturb NGC 4330. The 
nearest large galaxy (Br = 12.5), the spiral NGC 4313, has a 
projected distance of 130 kpc and a line-of-sight velocity dif- 
ference of 126 km s^'. It could be close enough in space and 
velocity to be physically associated with NGC 4330, possibly 
causing tidal interactions, though the old stellar disk of NGC 
4313 appears regular and undisturbed. Since NGC 4330 is in 
a cluster environment, we also cannot rule out the possibility 
that it has been affected by fast gravitational encounters with 
other galaxies which are now located far away. However, such 
encounters could not be responsible for the highly disturbed 
ISM of NGC 4330, since its stellar disk appears undisturbed. 

3. OBSERVATIONS 
3.1. Optical Imaging 

Exposures of NGC 4330 were taken in Harris BVR and 
narrow-band Ha + [Nil] (W16) filters at the 3.5 m WIYN 
telescope at KPNO on February 3, 2006. The 4k x 4k OPTIC 
CCD was used, which has a plate scale of 0"14 pixel^', giv- 
ing a field of view of about 9'.3, or 45 kpc at Virgo distances. 
The narrowband filter had a bandwidth of 72 A about the cen- 
tral wavelength of 6618 A. Integration time totaled 9 minutes 
each for R and V, and 15 minutes for B and Ha + [Nil], using 
three exposures per filter The seeing in R was 0"7, and the 
seeing in BV and Ha + [Nil] was 1"1. Standard reductions 
were performed using IRAFQ- A continuum-free Ha + [Nil] 
image was produced by subtracting the sky-subtracted, scaled 

' IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation. 
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Fig. 1 .— NGC 4330 HI contours (red) on a ROSAT X-ray (Bohringer et al. 
1994) grey scale and contour map (black) of the Virgo cluster The galaxy has 
been enlarged by a factor of 15. Note that the tail roughly points away from 
M87, indicating that the galaxy is probably falling in towards the cluster core. 
The aiTows indicate the possible range in direction of motion, as infen'ed 
from the HI tail angle, the location of the galaxy's radio deficit regio n, and 
the position angles of dust features (see also Figure|2]and Section|^7). 
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Fig. 2.— BVR-Ha image of NGC 4330 with radio deficit region (red), HI 
C-an'ay contours (white), with BVR and Ha (pink). Pixels with a ratio of 
observed-to-modeled 20 cm radio continuum flux density < 50% are said to 
have a radio deficit (see Murphy et al. 2009). 

R-band image from the narrow band image. We also include 
a 4-color image (Figure |3]l to illustrate the galaxy-wide dust 
distribution and how it relates to the older stellar disk and Ha. 
This was constructed by first convolving the R-band to the 
same seeing as our B, V, and Ha images. 

3.2. HI 21 cm Line and 20 cm Continuum Observations 

HI 21 cm line and 20 cm radio continuum observations 
were taken as part of the VIVA (VLA Imaging survey of 
Virgo galaxies in Atomic gas) survey (for futher details see 
the VIVA atlas paper, Chung et al. 2009), using the C-short 
and D-array configurations. Observations were taken on Au- 
gust 8, 2005, and on December 1, 2005, respectively. The 
observation time was 8 hours in C-array and 2 hours in D- 
array, with a bandwidth of 3. 125 MHz centered at a heliocen- 
tric velocity of 1565 km/s. The velocity resolution of both 
cubes is 10.4 km/s. The synthesized beam size FWHM for 
the C-short array data is 16.5 x 15.9 and the FWHM for the 
D-array data is 72.2 x 48.2 Each data set has its own advan- 
tages: the C-short array is less sensitive to extended emission, 
but its superior angular resolution enables more meaningful 
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Fig . 3 . — BVR-Hct image of NGC 4330 with insets of interesting dust features. Eacli inset lias liad its contrast adjusted to best liigliliglit the dust features, and 
the approximate position angle of the linear dust features is indicated with a black line. In a) and b), we detect dust features which are elongated roughly parallel 
to the HI tail and perhaps the projected ICM wind direction. There are also several bright HII regions, which may be in the process of being stripped of obscuring 
dust. 



comparisons between the spatial distribution of HI and the 
better-resolved maps at other wavelengths. In order to com- 
bine the angular resolution of the C-array and the sensitivity 
of the D-array, we have constructed a C-i-D-array cube with 
a FWHM of 26.4 x 24.0 ". In this paper, we use the C-array 
data to compare the spatial distribution of HI with other wave- 
lengths. We use the C-i-D array data for calculating the total 
HI flux and mass. The data were calibrated and continuum 
subtracted using the NRAO Astronomical Image Processing 
System (AIPS). For further details on radio image processing 
see Chung et al. (2009). 

3.3. GALEX UV Photometry 

We imaged NGC 4330 with GALEX (Martin et al. 2005) 
with the NUV broadband (1750-2800 A) and the FUV broad- 
band (1350-1705 A) filters on May 13, 2007, with exposures 
reaching a depth of 1 .7 ksec in both bands (GI Program Num- 
ber 79, Cycle 1). The FWHM is 5.3" for the NUV filter and 
4.2" for the FUV filter, and the pixel size is 1.5 arcsec. Data 
were processed using the standard GALEX pipeline (Morris- 
sey et al. 2007). 

4. RESULTS 

The undisturbed outer stellar disk, traced in R-band (Fig- 
ures |4}; and|5]l, indicates that the galaxy has not been signif- 
icantly affected by gravitational processes, leaving hydrody- 
namic processes as the likely cause of NGC 4330's truncated 
star-forming disk and disturbed, asymmetric HI distribution. 
Ample evidence indicates ongoing ICM-ISM stripping, which 
is expected to remove a galaxy's ISM from the outside in and 
produce a one-sided tail of extraplanar gas while leaving the 
old stellar disk intact. Ram pressure may also initiate a burst 
of star formation along the ICM-ISM interaction boundary 
and/or in the stripped gas. Below, we describe the evidence 
that various components of the galaxy's ISM are experiencing 
ram pressure stripping. 

4.1. The Old Stellar Disk 

The inner disk of the galaxy shows modest asymmetries in 
R-band light due to asymmetries in star formation and dust 
extinction, but the outer stellar disk beyond the dust and star 
formation (^75" from the galactic center) is quite symmet- 
ric and regular The R-band disk remains symmetric out to 
the deepest extent to which we can detect emission (^ 11.8 
kpc (2.5'), or 110% of Rjs)- The outer R-band isophotes are 
roughly concentric (FigureHJ, and there is no evidence of any 
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Fig. 4.— a) NGC 4330 Ha + [Nil] greyscale image of NGC 4330- ongoing 
star formation is largely confined to the inner disk, with a small amount in 
the tail; b) NUV greyscale - recent star formation is slightly more radially 
extended than HI; c) HI C+D-an'ay contour on R greyscale - the older stellar 
disk is undisturbed, but HI has been removed from the galaxy's outer disk 
(contours are 2, 4, 8.... 256 x lO'^cm^^) 

high- or low-luminosity bridges, warps, or other features as- 
sociated with gravitational interactions. The rotational period 
at this radius is ^ 10^ years, so if the outer stellar disk had 
been significantly perturbed in the past few Gyr, the evidence 
would not yet have been wiped out by the galaxy's differen- 
tial rotation. Although there is a low-luminosity tail of young 
star clusters roughly coincident with the HI tail (see the deep 
B-band image in Fig. |6ll, there is no smooth tail of older stars 
which, if present, would be a signature of gravitational inter- 
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action. 

4.2. Dust Extinction Morphology 

We have used the R-band images for analysis of the 
galaxy's overall dust distribution due to its higher resolution 
in our data (0"7 seeing in R vs. 1"! in BV), although we 
also use B and V bands for our analysis of the morphology 
of individual features. On the southeastern side of the disk, a 
spatially large 10" or 800 pc in its longest dimension) and 
significantly obscuring dust plume is elongated almost per- 
pendicular to the galaxy's photometric major axis (Figure|6]l. 
This cloud seems to mark the truncation boundary of the dust 
extinction on the northeastern side of the galactic disk. Be- 
yond it, no extinction is visible in the BVR bands, even in the 
B-R image (Figurel?) which emphasizes the galaxy's dust dis- 
tribution while eliminating the gradient caused by the disk's 
radial decline in luminosity. We find that at roughly the same 
radius as the dust plume referred to in Figure|6]as the "upturn 
dust cloud," UV and Ha emission are also truncated, curv- 
ing away from the disk midplane and towards the large dust 
cloud, i n a feature we call the "upturn" (discussed further in 
Sectionl43]l. 

On the northeast (leading) side of the galaxy, there are at 
least two smaller-scale dust clouds with elongated morpholo- 
gies (Figure|3^, b). One end of each of these clouds is located 
upwind (northeast) of the galaxy's continuous dust distribu- 
tion, and they extend downwind (southwest) to the galaxy's 
main dust lane. The cloud in Figure [3^ is ^6" (480 pc) long 
before it becomes indistinguishable from the galaxy's main 
dust lane, and its maximum width is 1.5" (120 pc). The cloud 
in Figure[3]5 is ^7" (560 pc) long by 2.5" (200 pc). These may 
be dense clouds which decoupled from the stripp ed lo wer- 
density ISM and are now being ablated (see Section l477] i. 

The dust distribution on the trailing (southwest) side of the 
galaxy at the base of the tail also curves away from and out 
of the disk. However, it curves much more gradually than 
the sharp cutoff we observe at the "upturn". Dust is visible 
out to about 50% of R25 on the leading (northeast) side, and 
there is some dust extinction out to approximately 90% of Ris 
on the opposite, trailing (southwest) side, although it is offset 
to the southeast of the major axis beyond about 50% of /?25- 
No extended dust structures are visible northeast of the major 
axis. 

4.3. HI Distribution 

The HI distribution in NGC 4330 is radially truncated 
within the optical disk along the major axis and visibly asym- 
metric, with extraplanar gas and a tail on the southwestern 
side of the optical disk, and little extraplanar HI northeast of 
the major axis (Figure]?] bottom panel). The HI distribution 
is significantly truncated on the leading (northeast) side of the 
optical disk (Figure [TsTi, extending to 55% of Ros along the 
major axis (Figure [SJ. On the trailing (southwest) side, the 
HI distribution extends to 95% of R25 along the optical ma- 
jor axis, although the actual tail, which is located southwest 
of the major axis, extends further In Figure |2] we define the 
main body and tail of the galaxy, as observed in HI. 

There is a large asymmetry about the major axis in the main 
body component of HI (region defined in Figure |2]l. The HI 
distribution peaks downwind of the optical major axis, and 
there is a large excess of emission downwind (southeast) of 
the galaxy (Figure|9|l. There is little doubt that most of this HI 
is extraplanar since it is located far outside the disk (Figure 



|4j;), and the HI distribution's asymmetry about the major axis 
allows us to estimate quantitatively the amount of extraplanar 
gas. We find that 22% (1.4x lO^M©) of the gas in the "main 
disk" component of NGC 4330 is actually extraplanar (Figure 
|9ll. 10% of its total HI located in the tail, so a full 32% of the 
galaxy's total HI emission is extraplanar. 

4.4. HI Kinematics 

NGC 4330 is redshifted by 469 km/s with respect to the 
Virgo cluster's systemic velocity, assuming a systemic ve- 
locity of 1100 km s^'. This is not an extreme value com- 
pared to other Virgo galaxies, and given that the galaxy has 
an HI tail pointing roughly away from M87 (Figure [Til, it is 
likely that there is a significant plane-of-sky component in the 
^northeast direction. Thus, the galaxy's plane-of-sky motion 
probably has a component towards the cluster center If NGC 
4330 has a typical 3D Virgo Cluster velocity of -1000 - 1500 
km s^^ then its plane-of-sky component is ^^700 - 1300 km 
s-i. 

In Figure [TO] we indicate the regions used in our compar- 
ison of the kinematics of planar and extraplanar HI. The HI 
along the "major axis" slice (Figure (TT] top panel) has rela- 
tively symmetric and normal kinematics, which suggests that 
the kinematics of most of the gas in the disk has not been 
strongly affected by ram pressure. The HI along the "down- 
wind" slice is located in the region of the g alaxy where we de- 
tect significant excess HI (see Section |43] l. and its kinematics 
show signs of disturbance. 

We find that the extraplanar gas on the approaching side of 
the disk (Figure (TT] bottom panel) shows a modest (~ 10 - 
20 km/s) systematic offset towards less extreme galactocen- 
tric velocities compared to the gas along the major axis. This 
differs from some galaxies undergoing active stripping (e.g. 
NGC 4522, Kenney et al. 2004, Vollmer et al. 2006), where 
the extraplanar gas has a pronounced overall displacement 
towards the mean cluster velocity. Since NGC 4330 is red- 
shifted with respect to the cluster mean, we might expect the 
stripped gas on the approaching (southwest) side to be accel- 
erated towards us with respect to the galaxy's rotation as it is 
swept towards us by the ICM along the galaxy's direction of 
motion. It is possible that the ICM wind has not affected the 
stripped gas close to the disk enough to accelerate it to more 
extreme velocities and the ccgas has lost angular momentum 
during the stripping process, lowering its galactocentric ve- 
locity. 

Figure [TT] (bottom panel) shows that the extraplanar gas on 
the leading side (northeast of the galactic center) shows lit- 
tle systematic velocity displacement. However, it does show 
smaller-scale disordered motions. We have identified several 
kinematic features of the extraplanar HI which are robust even 
to high flux cutoffs, labeled A-D in Figures [10] and [TT] bot- 
tom panel. In Figure [TT] bottom panel, feature A indicates 
that some gas is redshifted away from the cluster velocity at 
about 30" from the galaxy center, while feature B shows gas 
blueshifted by —10 km/s towards the cluster velocity at a very 
similar radius. Feature C also shows gas blueshifted towards 
the cluster velocity. These disordered motions suggest turbu- 
lence in the stripped extraplanar gas on scales of —20" (1.6 
kpc) and —25 km/s. 

A map of the HI velocity dispersion is shown in Figure [T2| 
The peak of the velocity dispersion is located near the galac- 
tic center, as is true for most normal spirals. The velocity 
dispersion in the tail of NGC 4330 is relatively low, peaking 
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Fig. 5. — NGC 4330 R contours (convolved with a Gaussian kernel to 3" resolution) on Ha greyscale. The outer R-band contours are roughly elliptical and 
regular (except for point sources). Each contour level represents a factor of 2 increase in flux from the previous contour, and the outermost contour is at 1.1 R25 
(25.6 mag arcsec^^). 
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Fig. 6. — NGC 4330 deep B-band image. Note the star clusters in the tail region and the upturn dust cloud. 
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NGC 4330 B - R image, where darker areas indicate relatively strong B-band emission and light areas indicate relatively strong R. 
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Fig. 8. — Normalized flux along NGC 4330's major axis in HI (C-airay), 
R, Ha, NUV, and FUV. Images were convolved to 16"resolution to match the 
HI. The flux is measured in a 15" strip along the major axis, and binned in 5" 
increments, then normalized. Error bars result from Poisson and sky noise. 
Compared to R, Ha is strongly radially truncated but symmetric about the 
galactic center. HI is slightly less symmetric, and more extended on the tail 
side. NUV and FUV are very asymmetric, indicating rapid evolution on 10 - 
100 Myr timescales on the tail side of the galaxy. 

at only '--^lO km/s. Except for the galactic center and the tail, 
the spectral line width is the same or higher in the extrapla- 
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Fig. 9. — Normalized flux along the minor axis of NGC 4330, measured in 
5" bins, where kpc indicates the major axis. HI emission peaks ^0.5 kpc 
downwind (SE) of the galactic center, R and Ha peak at the center, and FUV 
and NUV peak ~0.5 upwind (NW). 

nar gas than the major axis gas at a given radius. However, 
this effect is modest and may be explained by greater velocity 
gradients in the extraplanar gas, rather than indicating a true 
increase in local velocity dispersion. 
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Fig. 10. — HI velocity greyscale map of NGC 4330 overlaid with HI veloc- 
ity contours. The regions labeled "major axis" an d "d ownwind" correspond 
to the the position- velocity information in Figure [TTI The cross mark indi- 
cates the kinematic center, with a velocity of 1568 km/s. The contours are 
at intervals of 20 km/s. Key kinematic features A-D are labeled here and in 
Figure rm bottom panel. 



4.5. The Ha and UV Distributions: Star Formation Tracers 

There is a pronounced asymmetry along the major axis - 
on the upwind (northeast) side, UV emission is truncated at a 
significantly smaller radius than on the downwind (southwest) 
side (Figure [T3Tl. There is also a modest excess of flux at Ha 
and UV wavelengths downwind (southwest) of the galaxy, in- 
dicating star formation in the stripped material. The observed 
asymmetries are strong evidence that the galaxy is being more 
powerfully affected by the ICM wind on the leading (north- 
east) side than the trailing (southwest) side. 

The UV in the disk plane extends much further on the trail- 
ing side than the leading side (Figure[8]i, with a more contigu- 
ous tail structure than is present in Ha. On the leading side, 
NUV and FUV emission peak at 5 kpc and have dropped by 
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Fig. 1 1 . — Top: Spatial velocity diagram of NGC 4330's major axis, de- 
fined as +I- 7.5" about the galaxy's major axis. Contours in both panels are 1 
mJy/beam x 2, 4, 8, 16, 32. Bottom: Spatial velocity diagram of major axis 
and extraplanar material. The greyscale is the same as in the top panel, and 
the contours ai'e +/- 7.5" centered 20" downwind (southwest) of the major 
axis (see labels on Figur e|10| for beam locations). Key kinematic features 
A-D referenced in Figure flOl are labeled here. Contours are same as top. 
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Fig. 12. — Top: HI velocity dispersion greysc ale and contours for NGC 
4330. Labels refer to same regions as in figure [Tol Contours are 10, 15, 
20, 25, 30, 35, 40 km/s. Bottom: Same greyscale, HI isovelocity contours. 
Contours are 1560 km/s ±20, 40, 60... 

80% at 6 kpc. On the trailing side they drop more gradually, 
with some faint emission extending to or just beyond R25. On 
both sides of the galaxy, UV emission tapers off more gradu- 
ally than the Ha. Emission at all observed wavelengths other 



than UV is much more symmetric. Along the major axis, the 
Ha flux is truncated at the same radius on the leading (north- 
east) and trailing (southwest) sides of the center, near 50% of 
^25 (Figure[8]l. It terminates at the same radius on the leading 
side as the dust extinction and the HI emission (Figure [T3]l. 
However, on the trailing side, HI extends a bit further along 
the major axis than Ha, and the dust extinction does not have 
as distinct a boundary. The Ha emission is less extended than 
the UV both radially and in the z-direction. The fact that the 
UV distribution is much more asymmetric than the Ha and 
HI along the major axis suggests that the ISM distribution has 
recently changed. 

The Ha emission is closely confined to the disk plane, but 
the UV distribution is more extended both upwind and down- 
wind of the disk. (We discuss extraplanar star formation in 
discrete regions downwind of the galaxy in Section 14.101 but 
in this section we include both these regions and the smooth, 
continuous emission between 1-4 kpc from the major axis.) 
In Figure [141 we show the UV and Ha distributions relative 
to the major axis. Of the Ha emission outside of the tail, 
89% is located within -I-/-0.5 kpc of the major axis (between 
the two red lines in Figure [T4l i. The Ha distribution is the 
least extended of all the wavelengths we examine (Figure |9]l. 
Its brightness peaks near the galaxy's center rather than up- 
wind, like the UV, or downwind, like the HI. Excluding the 
tail, only 0.3% of the Ha is located more than 1 kpc away 
from the major axis (between the two blue lines in Figure 
[Hi, compared to 15% (18%) of the FUV (NUV). Emission 
far from the galaxy due to foreground or background sources 
has been masked out, and we count only pixels with a surface 
brightness greater than 3o above background, and count only 
emission from the main body, not the tail. 

Of the UV emission that is located 0.5 - 1 kpc from the 
disk plane, there is a significant excess of NUV and FUV 
flux downwind of the galaxy - 12% of the galaxy's main body 
FUV and NUV emission is 0.5 - 1 kpc to the downwind side, 
versus 3% (6%) of the FUV (NUV) upwind. The flux upwind 
is likely attributable to the radial extent of the UV within the 
disk, when viewed at a disk inclination of of 84°, and is prob- 
ably not extraplanar. The excess downwind FUV flux and the 
high FUV/NUV ratios of the extraplanar starlight indicate re- 
cent star formation in the stripped gas. 

As we discuss in Section 14.101 several discrete UV-bright 
regions are detected downwind of the galaxy at distances be- 
tween 1.5-9.5 kpc from the major axis, and most of their UV- 
determined stellar population ages are consistent with popu- 
lations too old for detectable Ha emission. However, it is also 
possible, given the ages of some of the regions, that we are de- 
tecting young stellar populations traced by UV but not Ha. It 
is well-established (e.g. ThiUcer et al. 2007) that many spiral 
galaxies have UV emission without associated Ha emission in 
their outskirts, and this may be due to the same physical con- 
ditions that are at work in the extraplanar UV regions in NGC 
4330. The outskirts of spiral galaxies that have UV emission 
but lack Ha are generally lower-density regions than the in- 
ner parts of the galaxies, and several possible explanations 
arise from that fact. The ionizing photons may leak out of 
the star-forming regions before creating HII (Oey & Kenni- 
cutt 1997). Alternatively, the IMF may be undersampled in 
low-density regions, leading to a shortage of the high-mass 
O and B stars necessary to create HII regions (e.g. Weidner 
& Kroupa 2006), or physical conditions may actually prevent 
massive stars from forming in low-density regions (Krumholz 
& McKee 2008). 
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Fig. 13. — Normalized flux along NGC 4330's major axis in HI, R, Ha, NUV, and FUV. Ha and HI emission on the leading side are truncated at almost the 
same radius as the dust extinction. 




Fig. 14.— NGC 4330 NUV, FUV, and Ha, with emission less than 30 
above background and emission unassociated with the galaxy masked out. 
By excluding the bright central region of the galaxy, we detect excess flux 
downwind in these wavelengths, a sign of recent and ongoing star formation 
in the stripped gas. 



4.6. The Upturn of ISM and Young Stars at the Leading Edge 

NGC 4330 has a remarkable feature at its leading edge 
which we refer to as the "upturn", shown in Figure[T5] at the 
northeastern Ha truncation boundary, the Ha emission veers 
abruptly out of the disk midplane at a ^ 35° angle. The Ha 
upturn has a vertical extent of '--^ 13" ( 1 kpc) from the optical 
major axis (Figure O. The large plume of dust described in 



Section|42]begins at the downwind edge of the Ha emission 
and extends about 0.5 kpc further downwind. We include a 
schematic of the upturn structure (Figure [TSI l to illustrate the 
relative positions of the different components. The FUV and 
NUV distributions gradually taper off in the radial direction 
beginning at the edge of the Ha. A hint of the upturn struc- 
ture is visible in the HI contours (Figure [TsT i despite the low 
HI spatial resolution. 

The upturn provides a window into the stripping process of 
the complex, multi-phase ISM. In Figure |7] made by dividing 
the B-band image by the the R-band image, we show that the 
B-to-R ratio is relatively large on the upwind (northeast) half 
of the upturn, and relatively small on the downwind (south- 
west) half. B is more susceptible to dust extinction than R, and 
is a better tracer of young stellar populations. The strongly 
blue color of the upwind (northeast) half of the upturn struc- 
ture (indicated in Figure[T6]as candidate star clusters) strongly 
suggests the presence of young stars in this area. Although 
this region contains some of the bluest B-to-R ratios in the 
galaxy, the total luminosity of this region is relatively small, 
and the color could be caused by a relatively small number 
of young stars. The redder areas in the B-to-R image are lo- 
cated on the downwind (southwest) side of the upturn, and are 
coincident with the large optical dust plume (shown in, for in- 
stance. Fig. |6]l, and the dividing line between high B-to-R 
and low B-to-R in the upturn may represent the boundary be- 
yond which the obscuring dust has been stripped away from 
the young star clusters. 

A key as-yet unanswered question about the stripping pro- 
cess is whether the denser parts of the ISM, such as the molec- 
ular clouds that could potentially host star formation regions, 
are stripped from the disk plane at the same time as the less- 
dense components such as HI. We detect no Ha emission or 
dust extinction in the disk at radii beyond the upturn, which 
strongly suggests that at most a very small fraction of the ISM 
is left behind in this part of the galaxy. A t the same time, the 
linear dust features described in Section 14.21 are located up- 
wind of the galaxy's main dust truncation boundary, and their 
morphology indicates that they may be in the process of being 
ablated by the ICM wind. The existence of these dust features 
suggests that although dense clouds can be left behind after 
the less-dense phases of the ISM are stripped, the remaining 
clouds are short-lived, which could be why we do not detect 
any in the disk beyond the upturn. 

We think that the Ha upturn contains stars which are form- 
ing in gas that has been pushed somewhat downwind of the 
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Fig. 15. — The northeast (leading) side of NGC 4330, including the "up- 
turn" region: a) HI C-an'ay contours on Ha - the Ha emission is truncated 
within the HI, and the Ha distribution at the galaxy's leading edge curves 
downwind with a hook-like moiphology; b) HI on NUV - the UV emission 
extends further radially than the HI and the Ha; c) HI on BVR-Ha - the R- 
band stellar disk extends well beyond the emission at other wavelengths, but 
the dust extinction is confined to much smaller radii. Note the large, sig- 
nificantly obscuring dust cloud which is located inside the third HI contour. 
Contour levels are 2, 4, 16 x lO" cm"^. 
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Fig. 16. — Schematic of the upturn region in NGC 4330. 



disk plane by the ICM wind. This would explain the up- 
turn's hook-like shape, as gas further from the galactic cen- 
ter has been displaced further downwind than gas closer to 
the center. We interpret the extended dust plume as ISM that 
has been moved even further downwind, perhaps because it is 
more susceptible to stripping - it could be more diffuse than 
the gas that was sufficiently dense to form the stars that now 
emit Ha in the upturn region. 



4.7. ICM Wind Angle 

The relationship between the true wind angle a, the angle 
between the ICM wind and the galaxy's disk plane, and the 
observed, projected angle 5, the difference between the posi- 
tion angle of the projected ICM wind direction and the posi- 
tion angle of the galaxy's major axis, is given by 



sina = sin5sin/sin(|) + cos/cos(|) 



(1) 



where (|) is the angle between the ICM wind direction and the 
line of sight and / is the galaxy disk's inclination relative to 
the line of sightQ In principle we can observe / and 5, and 
we can constrain (|) using the galaxy's measured line-of-sight 
velocity: 



tan(|) 



(2) 



where Vj^, is the wind component in the plane of the sky, and 
v/os is the wind component along the line of sight. In the case 
of a static ICM, these are equal and opposite to the galaxy's 
velocity in the plane of the sky and along the line of sight, 
respectively. 

In the case of an edge-on galaxy with an inclination of 90°, 
the second term in Equation 1 vanishes, and the equation re- 
duces to 



sin a = sin 5 sin (|) . 



(3) 



We identify three observational constraints on the projected 
ICM wind angle 5 in NGC 4330: the location of the local ra- 
dio deficit region, the angle of the HI tail, and the position 
angle of several smaller-scale, elongated dust features. The 
position angle of the large-scale tail will be affected by the 
stripping history over the past several hundred Myr, while the 
radio deficit region and the small-scale dust features proba- 
bly reflect more recent stripping conditions. When taken into 
account, NGC 4330's inclination of 84° can only potentially 
decrease a by ^^5°, which we will show is less than the error 
bars on our measurements. We therefore treat the galaxy as 
purely edge-on in the following analysis. 

The local radio continuum deficit region in NGC 4330 
(Murphy et al. 2009) is located approximately opposite the 
HI tail, to the north (upwind) of the galaxy's major axis (Fig- 
ure |2]l. Murphy et al. (2009) find that radio deficit regions 
are located near the leading edges of galaxies experiencing 
ram pressure. In the idealized case of a spherically symmet- 
ric ISM, the projected wind angle should be equal to the an- 
gle Qro, the angle between the galaxy's major axis and a line 
drawn from the galaxy center through the centroid of the ra- 
dio deficit region. For a disk galaxy, the ram pressure is more 
easily able to affect the low density gas (and weaker magnetic 
fields) in the halo than the high density gas (and stronger mag- 
netic fields) in the disk (e.g. Vollmer et al. 2009). Thus it is 
easier to create a radio deficit in the halo than the disk, and as 
a result the angle between the galaxy center and the centroid 
of the radio deficit region is offset toward the minor axis (the 
halo-dominated region) from the projected wind angle 5. De- 
termining this offset quantitatively would require modelling 

^ An analogous problem is presented in Kinney et al. (2000), who describe 
the relationship between the observed and projected orientation of the disk of 
a Seyfert galaxy and its jet. The geometry is nicely described in their Figure 
2 and Table 4. We substitute the ICM wind for the Seyfert jet, and define our 
angles the same way with the exception a, which we define as 90° — P 
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beyond the scope of this paper. For this paper, we give rough 
estimates. 

The radio deficit region curves around the disk on the NE 
(leading) side, extending all the way to the major axis. The 
curvature shows that in this region, a radio deficit has been 
created in material with higher density and stronger magnetic 
fields than those parts of the galaxy further from the major 
axis where most of the radio deficit region is located. Us- 
ing our definition in the previous paragraph, Qrd = 60 ±10°. 
Since 5 < Qrd, as explained above, the projected angle 5 is 
< 70°. The curvature of the deficit region means that pro- 
jected angles down to ^30° with respect to the disk plane are 
possible, so the radio deficit region's morphology suggests a 
projected wind angle between 30° - 70°. 

The position angle of the HI tail is another constraint on the 
ICM wind angle. In the case of a static ICM, the centroid 
of the trail of stripped gas will correspond to the galaxy's 
orbit through the cluster (e.g., Roediger & Bruggen 2007). 
In models, the relationship between the projected wind an- 
gle and the tail angle is influenced by many factors, including 
the length of the tail being considered, the average density 
of the tail gas, the evolutionary stage of stripping, whether 
the model uses hydrodynamics or sticky-particle n-body sim- 
ulations, and which parts of the tail are used to measure the 
angle. In certain cases, the galaxy's tail is fairly well-aligned 
with the ICM wind direction. Model snapshots in Roediger 
& Bruggen (2006) show that in early-stage stripping (^200 
Myr after the start of pressure), the high-density gas near the 
trailing side of the disk forms a tail ^-^20 kpc long whose po- 
sition angle is within ^15° of the ICM wind angle. The two 
sides of the galaxy are only asymmetric when the ICM wind 
is inclined from face-on, so the existence of distinct leading 
and trailing sides in NGC 4330 is consistent with a wind angle 
that is inclined from face-on. No systematic studies of the ef- 
fect of ICM wind angle on asymmetries in stripped gas exist, 
but from snapshots in Roediger & Bruggen (2006) and Roedi- 
ger et al. (2006), the asymmetry is apparent in galaxies with 
a wind angle of 60° (with 90° being face-on stripping). In 
other papers they treat pure face-on stripping, but a 60° angle 
is the largest non-face-on angle they consider, so we do not 
know the minimum offset from a face-on wind that produces 
an identifiably asymmetric HI distribution. 

As stated in Section|4j9] the HI tail is oriented at a 45 ±5° 
angle to the disk plane. It takes time for stripped ISM to 
move downwind, so we would expect the position angle of 
the denser tail gas close to the disk to increase over time, 
and ultim ately approach the ICM wind direction (e.g. Sec- 
tion |4|9]l. For a constant ICM wind angle, it seems unlikely 
that the HI tail near the disk would form a larger angle with 
the disk plane than the ICM wind does - to first order, the 
remaining gas disk should shield the recently-stripped mate- 
rial just downwind of it. Simulations consistently show that 
the tail's position angle does not exceed the ICM wind an- 
gle in the absence of large-scale turbulence or fallback of the 
stripped material. The importance of these complicating ef- 
fects is not yet well-understood, but to first order, we interpret 
the tail's position angle as an approximate minimum projected 
ICM wind angle. As previously stated, we do not know the 
precise minimum inclination with respect to the ICM wind 
that produces asymmetric stripping, so we estimate an upper 
limit on the projected ICM wind angle of ^85°. 

Assuming that the local gas flow direction within the galaxy 
is the same as the mean ICM wind direction, the position an- 
gle of small, dense clouds being ablated by the ICM wind 



are likely an indication of the projected wind direction. The 
small, elongated clouds in Figures [3^ and I3J), described in 
Section l4~2l have position angles of 35° ± 15°. There could 
be large systematic uncertainties in this method of determin- 
ing the projected wind angle, since the local gas flow within 
the galaxy may differ from the mean ICM wind direction, but 
the position angles of these dust clouds are consistent with the 
wind angle determined by the two other methods. 

In summary, the morphology of the radio deficit region in- 
dicates a projected wind angle of 5 = 30° - 70°, the HI tail 
is consistent with a projected wind angle of 40° - 85°, and 
the small dust features constrain the local flow direction at 
their location in the disk to 35° ± 15°. If all three measure- 
ments are considered, the projected wind angle is constrained 
to 40° — 50°. Given the unknown systematic uncertainties as- 
sociated with measuring the wind angle from the dust clouds, 
we prefer to constrain it based on only the HI tail and the radio 
deficit region, giving a projected wind angle 5 = 40° — 70°. 

The true wind angle a is a function of both 5 and (|), the an- 
gle between the true ICM wind direction and our line of sight. 
We can constrain (|) using Equation 2 and estimates for the 
wind velocity components v/o, and Vsky. The galaxy's line-of- 
sight velocity with respect to the Virgo Cluster's systematic 
velocity is H-469 km s^', so assuming a static ICM, the ICM 
wind velocity along the line of sight is y/„5 = -469 km s^' 
(towards us). Since the tail is elongated in the plane of the 
sky, the ICM wind also clearly has a component in that di- 
rection, although the magnitude of the plane-of-sky velocity 
Vsky- is hard to determine. Based on the morphology of the 
tail, Vsky > v/oi, so we adopt 469 km s^^ as the minimum rea- 
sonable value. No Virgo galaxies have known total orbital 
velocities > 2000 km s^^, so we adopt this as a maximum. 
Using Equation 3, we find that —75° < (|) < —45°, where the 
negative sign indicates that the wind angle is projected in the 
far hemisphere of the galaxy (due to the negative line-of-sight 
velocity). 

We can now constrain a using Equation 1 and our derived 
estimates of 5 and (|). Since sin(|) will never be >1, a < 5, 
meaning that the 3D angle between the disk and the wind will 
never be greater than the projected angle. Thus, a < 70°. A 
lower limit on a can be calculated using the maximum (|) and 
minimum 5 values, with the result that 30° < |a| < 7013. 

4.8. The Gas-Stripped Outer Disk: Estimation of Stripping 
Rate from UV-Optical Colors 

The leading side of NGC 4330 has a gas-stripped outer 
disk between 5-8 kpc from the galactic center where Ha 
and HI emission and dust extinction are faint or undetected, 
but where we can still detect both FUV and NUV emission. 
Ha emission is a tracer of very recent and ongoing (within 
^ 3 X 10*' years) star formation, and UV light traces young 
(NUV traces ^10*^ years old or less, FUV 10^ years or 
less) stellar populations (e.g. Bruzual & Chariot 2005), so 
the 5-8 kpc annular region which lacks Ha emission but 
emits strongly in the UV has had star formation quenched 
sometime within the past ^0.5 GyrQ Assuming that star 
formation was occurring at the time of stripping, using the 

^ The sign of a is negative in this case because, consistent with the negative 
value of (|), the wind first encounters the galaxy in the far hemisphere. 

As discussed in Section l431 UV emission without accompanying Ha is 
sometimes detected in the outer disks of spiral galaxies, and does not defi- 
nitely rule out ongoing star formation. However, the area of NGC 4330 near 
the Ha truncation radius is well within the disk, where we would not ex- 
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Fig. 17.— log(FUV/NUV) image of NGC 4330, where the input images 
have been convolved to 9" resolution to minimize noise. Identifiable back- 
ground sources have been masked. Boxes indicate the apertures used to mea- 
sure flux from to determine the relative UV colors on the leading (NE) and 
trailing sides (SW). The apertures are located in the stripped outer disk, be- 
yond most of the galaxy's dust and ISM. The stellar populations on the lead- 
ing side have a lower log(FUV/NUV) value, indicating that the stellar popu- 
lations there are older than the ones on the trailing side. 

FUV/NUV/optical flux ratios to measure how long it has been 
since the youngest stars formed in a stripped region will pro- 
vide an estimate of how much time has elapsed since the area 
was stripped and star formation was extinguished. Although 
UV is very prone to dust extinction, optical images of NGC 
4330 show little or no dust extinction on the leading side of 
the disk beyond the stripping boundary at 5 kpc. The trailing 
side shows some dust extinction along the major axis, which 
can be expected to redden the observed colors. By comparing 
the observed UV colors with stellar population models, we 
can potentially learn how long it has taken to strip the leading 
edge of the galaxy to its current radius and if the leading and 
trailing sides of the galaxy have different stripping timescales. 

We measured log(FUV/NUV) and log(NUV/r) in 15" high 
by 10" wide apertures along the galaxy's major axis from 5-8 
kpc (see Figure[T7]for aperture locations). We ran several dif- 
ferent Starburst99 (Leitherer et al. 1999) models in order to 
determine which type of stellar population best matches the 
UV-optical colors we observe. The single stellar population 
model simulates a stellar population that formed all at once 
in a burst, and is described by the time since that burst. The 
"quench" model reflects a constant star formation rate, fol- 
lowed by an abrupt truncation of star formation activity, and 
is described by the time since star formation stopped. The 
"quench with burst" model is similar to the quench model, 
except that it has a starburst of a given strength at the time of 
truncation. In Figure[T8] we indicate a model in which 5% of 
the region's total stellar mass was formed in a burst. 

In the gas-stripped outer disk, there is a large gradient in the 
UV-optical colors on the leading side, but virtually no gradient 
on the trailing side (Figure [TSll. At 5 kpc the colors are very 
similar on the two sides of the galaxy, with fairly "blue" col- 
ors indicating significant recent star formation. However, the 
colors "redden" significantly between 5-8 kpc on the lead- 
ing side, by ^0.8 AB magnitudes in FUV-NUV and ^0.9 in 
NUV-r This color gradient indicates that star formation has 
been gradually quenched from the outside of the galaxy in- 
wards over the course of a few hundred Myr on the leading 
side, and the lack of a similar gradient on the trailing side in- 
dicates that star formation was all quenched recently at about 
the same time. This difference in star formation histories 
reflects different stripping histories on the two sides of the 

pect to find the low-density star-forming regions necessary for ongoing star 
formation that lacks Ha emission. 



galaxy. Given the similar colors, the trailing side may have 
been stripped around the time when the gas at 5 kpc on the 
leading side was stripped. 

Several effects complicate the interpretation of the ob- 
served colors, including dust extinction, contamination from 
the outer disk due to the composite stellar population along a 
given line of sight, asymmetrical quenching times, and rota- 
tion. These are discussed in detail in the Appendix. Here we 
summarize the main results. 

After correcting for the effects of the composite stellar pop- 
ulation along each line of sight, we find that our data remain 
inconsistent with a single stellar population on the leading 
side (as expected for a disk galaxy), but consistent with star 
formation that has been gradually quenched from the outside 
in over the course of 200 - 400 Myr The UV/optical color 
gradient on the leading side is roughly parallel to the model 
age tracks, strongly suggesting that there is a stellar popu- 
lation age gradient from 5-8 kpc. We are confident that a 
gradient in quenching times exists in the stripped region on 
the leading side. If the gradient were simply a result of the 
viewing geometry we have attempted to correct for, we would 
see a similar gradient on the trailing side, but points mea- 
sured at the same radii on the trailing side all have about the 
same quenching time (Figure [Tsll. Since there is some dust 
extinction on the trailing side, we calculate a modest redden- 
ing vector of Afuv=l, which corresponds to E(B-V)=0.035, 
using the wavelength-dependent dust extinction formulae of 
Calzetti et al. (2001). This gives a quenching time of ^50 
Myr ago for the trailing-side regions. 

Just beyond the Ha boundary on the leading side, at 5.0 kpc 
from the galactic center, the corrected colors are consistent 
with a quenched stellar population with a quenching time of 
<50 Myr. The modeled burst fraction and quenching time in- 
crease with increasing radius. The maximum quenching time, 
reached at 8.4 kpc, is 200 - 400 Myr ago. Even after correc- 
tion, the burst fractions associated with the outer disk points 
are much larger than seems realistic, since ^-^20% - 30% of the 
current population would have been created in the burst. We 
suspect that uncertainties in stellar population modeling are at 
work. 

The recent work of Conroy et al. (2009, 2010 a and b) 
explores sources of uncertainty in stellar population model- 
ing such as the slope of the IMF, poorly-understood phases of 
stellar evolution, the effects of metallicity distributions within 
a population, and dust distribution on large and small scales. 
Changes in the detailed treatment of dust extinction or the 
IMF slope can each cause differences of >0.1 magnitude in 
modeled (NUV - r) colors (Conroy et al. 2010a). Although 
Conroy et al. do not analyze the Starburst99 model specifi- 
cally, the uncertainties they mention apply to all stellar pop- 
ulation synthesis models. We therefore use stellar population 
synthesis models with an appreciation of the inherent uncer- 
tainties. 

4.9. The Tail 

NGC 4330 has a long HI tail extending southwest beyond 
the optical disk which contains about 10% (O.4xlO**M0) of 
the galaxy's total HI. The HI tail is 13 kpc long and 5 kpc 
wide (Chung et al. 2007), and is oriented at a ^45° angle to 
the major axifl The tail structure is particularly interesting 

We measure thi'ough the middle of the highest contour of tail emission 
and assume that the vertex of the angle is located where the ta il b egins to 
bend away from the disk (x=0" and y=30" in all panels of Figure II 9t 
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Fig. 18. — FUV-NUV vs. NUV-r plot of the outer disk of NGC 4330. The dashed lines indicate different stellar populations models, from left to right: a 
quenched model in which constant star formation is abruptly truncated at the specified time, a similar quenched model that experienced a 5% starburst at the time 
of truncation, and a single stellar population model. The solid black lines labeled 50 - 400 Myr indicate a constant age or quenching time across the models. The 
points show the data before our correction algorithm has been applied, and the boxes show the range of parameter space which each point may occupy after our 
correction algorithm has been applied. The reddening vector has been calculated for Afi]v=i- There is a pronounced gradient in stellar population quenching 
times in the gas-stripped outer disk on the leading side (circles), where the approximate stellar population quenching time increases by ~300 Myr from 5-8 kpc 
in radius, but there is no gradient on the trailing side (squares). 



because there is a spatial offset between the HI and UV/Ha 
tails (Figure l20b. suggesting recent evolution, and this is dis- 
cussed below. 

There is some ongoing star formation in the tail, but it 
is sparse, consisting of several small Hll regions extending 
southwest from the main, continuous Ha disk. The tail ac- 
counts for about 2% of the galaxy's Ha flux, much less than 
the 10% of the total HI flux coming from that region. There 
is a large gap in between the emission at the base of the tail 
region and Hll regions at the tip of the tail (Figure [T9]l, which 
are located 2.6 kpc from the major axis. The tail's terminal 
Hll region is located just outside '^1.1 R2i- No smooth stel- 
lar component of the tail is visible in R, although the dust 
distribution on the tail side of the galaxy curves downwind 
(southeast) and out of the disk, and star clusters are visible in 
optical broadband in the tail region. 

The UV emission in the tail is stronger and much more con- 
tinuous than the Ha (Figure [T9k. b), with 8% of the galaxy's 
FUV and NUV flux coming from the tail. The UV tail has 
no gaps between the disk and tail emission, in contrast to the 
large gap seen in the Ha. The UV tail is oriented at a ^35° 
angle to the major axis, which is shallower than the HI tail 
(Figure l20b. although they both begin to bend away from the 
disk at about the same radius. The HI tail is much longer 
than the UV/Ha tail, extending nearly 3 kpc further from the 
disk. The sharp drop in UV surface brightness at the down- 



wind edge of the tail in Figure [T9b suggests that this is a real 
difference, and not simply a reflection of our sensitivity limits 
in HI and UV. 

The centroid of the HI tail emission is located downwind 
(southeast) of the UV tail, suggesting that much of the ISM 
has moved from the location of the UV tail to the current lo- 
cation of the HI tail. The C-array (highest-resolution) HI map 
(Figure [T9b) shows two main peaks in the tail HI distribu- 
tion, which are displaced from th e pe aks in the Ha (Figure 
[T9b ) and UV emission. In Section|531 we estimate that it has 
taken the HI gas 10 - 300 Myr to move to its current position 
downwind of the UV-bright stars. 

The lack of widespread Ha emission indicates that there 
has been little star formation in the past ^5 Myr The strong 
FUV emission and "blue" FUV-NUV colors (see Section 4.5) 
reveal that, between ^10 and 200 Myr ago, there was sig- 
nificant star formation activity in the tail. We can infer that, 
in the intervening time, conditions in the tail region became 
less hospitable to star formation - either the existing gas was 
used up by star formation, or it was swept away during the 
stripping process. The presence of significant amounts of HI 
downwind (southeast) from the stellar tail supports the latter 
hypothesis. 

4.10. Extraplanar Star Formation Regions 
4.10.1. UV Properties 
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Fig. 19. — NGC 4330 tail and southwest outer disk, a) HI C-array (highest 
resolution) on Ha - the tail in Ha is sparse, consisting of a few discrete re- 
gions along the northern border of the HI emission. Contour levels are 2, 4, 
16 X lO'^cm^^; b) HI on NUV - the NUV emission peaks are coincident with 
those in Ha, but there is also fainter, continuous emission connecting them 
with the main body of the galaxy; c) HI on BVR-Ha - the dust extinction in 
the disk decreases significantly at the same radius as the Ha; d) FUV on Ha - 
the Ha peaks ai'e located roughly along the centroid of the FUV tail, suggest- 
ing that the location of star formation has not changed in the last lO' — 10** 
years. 

We detect a striking number of extraplanar UV-bright re- 
gions downwind (southeast) of the disk plane, indicating the 
presence of recent or ongoing star formation in the stripped 
gas. There are 9 regions, at distances up to 9 kpc from the 
disk, which are distinct from the main body and tail of the 
galaxy and which have surface brightnesses of > 3o above 
background in either FUV or NUV. We have masked flux from 
objects whose optical morphologies reveal them to be stars or 
background galaxies. Although many of the extraplanar UV- 
bright regions are faint, they are all on one side of the galaxy. 
All but two of the UV-bright extraplanar regions are within 
the VLA HI contours (Figure l22li. and one of the remaining 




Fig. 20. — NGC 4330 schematic showing the distributions of R (red), FUV 
(green). Ha (purple), HI (blue), and 20 cm radio continuum (contours). The 
radio continuum and HI tails are displaced downwind (SE) of the UV and 
Ha tails, indicating that the ISM in this area has been pushed downwind over 
time. 

regions is close to the outermost detected HI. We measured 
the flux of each UV region using a circular aperture as shown 
in Fig. |2T|and give the results in Table 1 . In total these regions 
account for 5% of the galaxy's FUV emission, and 4% of its 
NUV. We detect Ha emission downwind of the galaxy from 
only 2 of the 9 regions, which together comprise only 0.3% of 
the galaxy's Ha flux. 
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Fig. 21. — Composite UV image of NGC 4330 where NUV is yellow and 
FUV is blue. We observe 9 extraplanar peaks of UV emission downwind 
of NGC 4330 which are not coincident with foreground stars or background 
galaxies in the optical images and presumably originate from recent star for- 
mation in gas stripped from the disk. The circular apertures used to measure 
the regions' photometry are shown. 

Most of the extraplanar UV-emitting regions are in two 
groups: one is a few kpc downwind (southeast) of the lead- 
ing edge, extending south beyond the Ha upturn (Figure |43] |. 
and the other is located downwind of the galaxy between the 
minor axis and the base of the tail. While the Ha upturn is 
a continuous curve of Hll regions that extends a modest dis- 
tance of about 0.8 kpc from the major axis, the UV beyond 
the edge of the upturn extends in a large, arc-shaped com- 
plex of regions that reaches ^7 kpc from the main body of 
the galaxy. These FUV-bright features may have formed from 
gas that was originally in the outer disk beyond a radius of 
~ 6 kpc (75"), a region which is now bright in NUV but faint 
in FUV. 

4. 10.2. Star Formation Histories of Extraplanar Regions 

The FUV/NUV ratios of these regions can be compared to 
stellar population models in order to estimate the amount of 
time that has elapsed since star formation stopped, as well 
as the stellar masses. We used Starburst99 stellar population 
models (Leitherer et al. 1999). We assume a scenario in which 
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star formation occurred in a short burst in pockets of stripped 
extraplanar gas, which is best described by a single stellar 
population (SSP) model. For the regions closest to the galactic 
disk (4,6,7, and 9), we have removed the contribution of the 
smooth disk by subtracting the flux in an aperture of equal 
size located opposite the major axis from the region itself. 

The FUV/NUV ratios and corresponding ages from the sin- 
gle stellar population model are plotted in Figure|23] Most of 
the regions are younger than 100 Myr Region 7 is the oldest 
at 350 Myr, and the second-oldest is Region 4, with an age 
of ~250 Myr. These are the two most massive regions, and 
it makes sense that we only detect old regions with relatively 
high masses - UV emission becomes relatively faint in older 
stellar populations, and there are likely other old regions with 
lower masses that are below our detection threshold. All of 
the other extraplanar regions have ages between '^lOO Myr 
at the oldest, and <10 Myr at the youngest. Regions 3 and 
5 are detected in FUV but lack significant NUV emission, so 
assuming an upper limit of 3o above background, we derive 
upper limits on their ages of ^10 Myr. The total stellar mass 
of the regions is 5 x IO^Mq, a tiny fraction of the galaxy's total 
stellar mass. 




Fig. 22. — Continuum-subtracted Ha with circles over distinct UV peaks, 
with HI C+D-array contours. Although UV emission is detected at significant 
(3a above background) levels in the Ha-bright tail. Ha is not detected at 
significant levels in any of the UV-bright regions away from the main body 
of the galaxy, with the exception of regions 4 and 8. HI contours are 2, 4, 
8,...xI0'''cm-2. 

The oldest of the UV-bright regions, which is one of the 
closest to the disk plane, is ^350 Myr old, and the quenching 
time of the ou ter d isk at the galaxy's leading edge is 200 - 400 
Myr (Section [4.8b . The removal of gas from the outer disk 
and the existence of relatively dense extraplanar gas are both 
consistent with strong stripping beginning around this time. 

4.10.3. Ha Properties 

We detect Ha emission from regions 4 and 8, to a sig- 
nificance of 4o and 8a above background, assuming point 
source morphologies. The detection of Ha emission from 
these two regions may seem somewhat puzzling, since they 
are not among the youngest regions according to their UV 
colors. While Ha emission can persist for up to 10 Myr in a 
single-burst stellar population, our UV-based age estimates of 
both of these regions are much older, 240 Myr ±10 Myr for 
Region 4 and 90 Myr ± 15 Myr for Region 8. 

We think the most likely explanation is that the regions 
with ongoing star formation may also contain dust that was 
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Fig. 23.— log(FUV/NUV) for the extraplanar regions in NGC 4330 indi- 
cated in Figures [2l1 and l22l plotted vs. distance to the galaxy's major axis. 
Regions 3 and 5 lack significant NUV emission, so their ages are lower lim- 
its. Regions 4 and 8 have associated Ha emission, and are probably not 
well-described by the single-burst stellar population model that provides the 
ages given on the right-side axis. 

Stripped from the galaxy along with the gas that is now form- 
ing stars. The star-forming region emits enough ionizing flux 
to be detected in Ha, but dust extinction reddens the UV color 
by preferentially absorbing FUV rather than NUV photons (as 
reflected in the wavelength-dependent extinction formulae of 
Calzetti 2001). An uncorrected FUV/NUV ratio wiU yield an 
age that is older than the region's actual age. 

The SSP models predict an Ha flux < 1 x 10"^^ erg s"^ 
cm^^ for the brightest region (Region 8) based on the UV- 
derived masses and ages, which is fainter than our 3o detec- 
tion limit for a point source. We do not detect Ha in Regions 
3 and 5, whose UV colors suggest that they are young enough 
to have lingering Ha emission from a recent episode of star 
formation, but their predicted fluxes are more than an order of 
magnitude below our detection limit. 

There are also other scenarios that could also explain UV 
emission without accompanying Ha in low-density or low- 
metallicity regions (e.g. Lee et al. 2009), so it may not neces- 
sarily be a sensitivity issue. As we discussed in Section |431 
UV emission without associated Ha is often detected in the 
outskirts of spiral galaxies (e.g. ThiUcer et al. 2005, 2007), 
which may indicate ongoing star formation with a very low 
surface density and/or a shortage of the O and B stars which 
create Ha emission. It is also possible that ram pressure strip- 
ping affects the star formation process to such a degree that 
standard assumptions about the IMF or the formation of HII 
regions are no longer valid. 

5. DISCUSSION 

5.1. The Evolutionary State of NGC 4330: Comparison With 
Other Virgo Spirals 

Of the seven Virgo spirals with long HI tails detailed in 
Chung et al. 2007, NGC 4330 displays the best evidence for 
strong, active ram pressure stripping. All of the galaxies in 
the sample have tails pointing roughly away from the cluster 
core, which is evidence that they are moving on significantly 
radial orbits towards the cluster center and are likely making 
their first pass towards the cluster core. Since they seem to 
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Extraplanar UV Region Properties 



Region 


FUV Flux 


FUV a 


NUV Flux 


NUV a 


Ha Flux 


Stellar Mass (M,.,) 


1 


2.80e-17 


3.4 


1.58e-17 


3.8 


<4.1e-16 


1.0e5 


2 


5.43e-17 


4.4 


2.97e-17 


4.7 


<4.1e-16 


8.2e4 


3 


2.50e-17 


3.5 


<1.08e-17 


<3.0 


<4.1e-16 


2.1e3 


4 


2.50e-16 


8.8 


1.70e-17 


12. 


8.58e-16(4a) 


2.1e6 


5 


1.69e-17 


3.5 


<7.48e-18 


<3.0 


<4.1e-16 


5.9e3 


6 


3.64e-17 


6.6 


I.99e-17 


9.1 


<4.1e-16 


3.9e5 


7 


5.49e-17 


7.3 


5.59e-17 


14. 


<4.1e-16 


1.6e6 


8 


1.13e-16 


10. 


5.88e-17 


10. 


2.60e-15 (80) 


1.6e5 


9 


6.21e-17 


11. 


4.35e-17 


14. 


<4.1e-16 


l.le5 



TABLE 1 

UV FLUXES ARE IN UNITS OF ERG S"' CM^'A"'. Ha FLUXES ARE IN UNITS OF ERG S"' CM"^. 



be falling inwards, they presumably have yet to reach peak 
ram pressure. The seven tail galaxies are all located near the 
cluster's virial radius, but NGC 4330 is the closest, at least in 
projection, to the cluster core, and may be encountering the 
densest ICM. NGC 4330 has the most sti'ongly truncated HI 
distribution and is the most Hl-deficient of the Chung et al. 
sample galaxies. 

The sample of Crowl et al. (2008) contains several Virgo 
spirals at more advanced stages of stripping than the Chung 
sample. Crowl et al. selected galaxies with star formation, 
as traced by Ha emission, truncated well within the stellar 
disk. The galaxies in their sample which are most likely to 
be analogous to NGC 4330 are the ones whose Ha distribu- 
tions are characterized as "truncated/normal" (Koopmann & 
Kenney 2004b), with no star formation beyond a truncation 
radius but normal star formation within that radius. These 
galaxies also have stellar disks well-fit by elliptical isophotes, 
with no evidence of asymmetries that might have been caused 
by gravitational interactions. Crowl et al. use stellar popula- 
tion analysis to learn the star formation history just beyond 
the Ha truncation radius, indicating when that region was 
stripped. The truncated/normal galaxies with young ( 50 - 300 
Myr) outer disk stellar populations include NGC 4522 (50 - 
100 Myr) and NGC 4569 (-300 Myr). Both galaxies show 
significant extraplanar Ha and HI, indicating relatively recent 
disturbances, and the HI within their disks is highly truncated 
to -0.3 R25. 

The HI truncation radius of NGC 4330 is further from the 
nucleus (-0.6 R25) than those of NGC 4522 and NGC 4569, 
implying that it is at the least advanced stripping stage of the 
three. NGC 4330 is less HI deficient than NGC 4522 and 
NGC 4569. NGC 4330 has -6 times less HI than expected 
for a galaxy of similar size in the field, and NGC 4522 has —8 
times less (Chung et al. 2009). NGC 4569 is much more HI 
deficient, with a factor of 30 less HI than expected. It also has 
a small HI tail pointing towards M87 (Vollmer et al. 2004), 
suggesting that it is post-peak ram pressure and therefore a 
later stage than NGC 4330 and NGC 4522. 

While NGC 4522 has been more strongly stripped of HI, it 
appears that stripping has been going on for a longer time in 
NGC 4330. The lai-ge color gradient in the outer disk of NGC 
4330 indicates that its stripping has occurred over a relatively 
long time, 200 - 400 Myr. In contrast, the relatively blue op- 
tical colors of the outer disk in NGC 4522 visible in SDSS 
images suggest that it was stripped during a recent, rapid, im- 
pulsive event. 

Truncated/normal galaxies with old outer disk stellar popu- 
lations (-400 - 500 Myr) in the Crowl sample include NGC 
4580 and IC 3392. These galaxies have highly truncated but 
symmetric HI with regular velocity fields that show no sign 



of recent disturbance. There is no evidence of extraplanar HI 
or star formation, and the galaxies are highly HI deficient. 
NGC 4580 has 32 times less HI than a similar undisturbed 
galaxy, and IC 3392 has 14 times less. These galaxies proba- 
bly underwent a similar process to the one currently affecting 
NGC 4330, but their higher HI deficiencies indicate that more 
of their HI has been stripped, and their lack of asymmetries 
or extraplanar material indicates that they experienced peak 
pressure at least several hundred Myr ago. 

5.2. The Evolutionary State of NGC 4330: The Orbit within 
the Virgo Cluster 

The evolutionary state of stripping depends on where the 
galaxy is along its orbit through the cluster Information on 
NGC 4330's orbit contributes to a consistent picture that it 
is at a pre-peak stage of strong stripping since it is moving 
inwards into the central region of the cluster, probably for the 
first time. 

The projected wind direction of PA = 17° 13° (Section 

14.7b measured counterclockwise from north (the number we 
give in Section |4j7]is the angle between the wind and the ma- 
jor axis) forms an angle of 45° — 75° with the line between 
the galaxy center and M87, where 0° is purely radial motion 
and 90° is purely tangential motion. If the ICM is static, the 
ICM wind direction indicates the galaxy's direction of mo- 
tion on the sky. Thus, the radial and tangential components 
of projected motion are similar. While there is a large tangen- 
tial component that will carry the galaxy to the north of M87, 
there is also a large component which is radially inward. 

We can estimate the galaxy's total motion through the clus- 
ter The galaxy has a line-of-sight velocity of 469 km s^ ' with 
respect to the V irgo cluster's mean velocity, and we estimate 
in Section Wn\ that the galaxy's plane-of-sky velocity is be- 
tween 500 km s^' and 2,000 km s^^ This gives a total cluster 
velocity between —700 - 2,100 km s~ '. B ased on the strip- 
ping history of the outer disk (Section 14.81 1. gas in the outer 
disk at radii greater than 8 kpc had been removed by —200 
- 400 Myr ago. Given the —55° projected angle between the 
galaxy's orbital path and the line to the cluster center, over 
a timescale of 200 - 400 Myr, the galaxy has moved 240 - 
490 kpc closer to the center. This means that strong stripping 
would have started when the galaxy was between 840 kpc - 
1.1 Mpc from the cluster center In the sample of galaxies 
with HI tails in Chung et al. (2007), the most distant in pro- 
jection from the cluster center is 0.9 Mpc. NGC 4330 could 
have been at this projected distance or closer when its outer 
disk was stripped if the stripping timescale is closer to 200 
Myr, but the timescale could be 400 Myr as long as its veloc- 
ity towards the cluster center has stayed under — 700 kms^^ 
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It is also possible that ICM substructure led strong stripping 
to start in NGC 4330 at larger cluster radii than the other tail 
galaxies - Virgo has significant substructure, and in a dynamic 
ICM, the ram pressure a galaxy experiences at a given cluster 
radius can vary by at least an order of magnitude (Tonnesen 
& Bryan 2008). Evidence for the dynamic ICM enhancing 
stripping has been found in another Virgo galaxy, NGC 4522 
(Kenney et al. 2004), which seems to be experiencing ram 
pressure an order of magnitude stronger than would be pre- 
dicted by static ICM models for its location in the cluster. 

5.3. The Evolutionary State of NGC 4330: Comparison With 
Simulations 

A comparison with simulations supports the idea that the 
ram pressure in NGC 4330 has increased significantly over the 
last -300 Myr. The ISM morphology of NGC 4330 is strik- 
ingly similar to the "medium ram pressure" model scenario 
of Roediger & Bruggen (2006) (see their Figures 5 and 7), 
with Pram 6.4 X 10"'^ erg cm"^\ t--200 Myr after the start 
of ram pressure, and an unprojected wind angle of 30° (their 
Figure 5) or 60° (their Figure 7) between the disk plane and 
the wind direction. These wind angles are roughly consistent 
with the --40° - 70° range we estimate for NGC 4330 (Section 
14.7b . The model is similar to the stripping history we propose 
in Section l4~8l which entails a significant increase in pressure 
200 - 400 Myr ago. The leading-edge side of their model also 
resembles the leading (northeast) side of NGC 4330, with the 
recently-stripped ISM near the leading side of the disk elon- 
gated almost perpendicular to the disk plane, much like the 
upturn as we detect it in dust extinction. 

It is clear that the stripping process has affected the lead- 
ing and trailing sides of the galaxy very differently - there 
is a clear stellar population age gradient on the leading side, 
but on the trailing side, there is no age gradient and the stel- 
lar populations are younger. The asymmetry observed in this 
galaxy is also consistent with a significant increase in ram 
pressure over the last —200-400 Myr. Both N-body (Vollmer 
et al. 2001) and hydrodynamic (Roediger & Bruggen 2006) 
models predict that as a galaxy is initially affected by ram 
pressure, the ISM distribution will become asymmetric be- 
cause the leading side is affected more strongly than the trail- 
ing side. Over time, asymmetries caused by stripping will 
even out as each part of the disk rotates through the region 
of highest ram pressure (Roediger & Bruggen 2006). Many 
details of the observed gas morphology of NGC 4330 are con- 
sistent with models with smoothly increasing ICM pressure, 
in which a galaxy falling into a cluster for the first time starts 
to experience strong ram pressure over the course of —200 
- 300 Myr (e.g., Roediger & Bruggen 2006, Vollmer et al. 
2001). 

All galaxies orbiting through clusters experience changes 
in ram pressure strength and direction. Steady changes are 
expected if a galaxy orbits through a smooth, static ICM, al- 
though changes may be more sudden in a dynamic, structured 
ICM. In NGC 4330, the UV disk asymmetry indicates a sig- 
nificant increase in ram pressure over the last —200-400 Myr. 
Over 400 Myr, NGC 4330 would likely move radially inward 
in the Virgo Cluster from —0.9 Mpc to —0.6 Mpc based on 
typical Virgo velocities. Over this time interval, the Virgo-like 
models (with a smooth, static ICM) of Roediger & Bruggen 
(2007) show increases in ram pressure by factors of 1 .3 - 2, to 
values of Pro,,, = 10^'^'^ — 10^^' '^erg cm^^ for galaxies un- 
dergoing pre-peak stripping. The wind angle can change by 



— 10° over 300 Myr for a galaxy which has not yet reached 
its closest approach to the cluster core, as is likely the case 
for NGC 4330. It is not yet clear whether the observed gas 
and stellar population distributions in NGC 4330 can be ex- 
plained by steady changes in ram pressure, or whether the 
sudden changes expected in a dynamic ICM are required. 

5.4. Tail Evolution 

We observe that the HI tail is offset to the southeast, down- 
wind of the UV-bright regions. Such offsets between gas and 
young stars can provide strong constraints on the recent strip- 
ping history of a galaxy. 

Spatial offsets between gas and recently-formed stars natu- 
rally occur with gas stripping (e.g., the models of Kronberger 
et al. 2008). Whereas the ISM continues to be accelerated by 
ram pressure, newly-formed stars decouple from the ISM, and 
retain the orbital energy and angular momentum from their 
birth. As the rest of the ISM continues to move downwind, it 
will leave a trail of newly-formed stars in its wake, resulting 
in a tail of young stars possibly oriented at a different angle 
from the disk plane than the HI tail. Small offsets between the 
newly-formed stars and HI are expected with small or nearly 
constant ram pressure, but large offsets, like those in NGC 
4330, may indicat e lar ge and changing ram pressure. As ex- 
plained in Section l531 the offset between the UV and HI tails 
may be consistent with with changes in the ICM wind angle 
and pressure that occur during a reasonable galaxy orbit. 

We can use stellar population ages of the UV-bright regions 
to place limits on the timescale of changes to the ISM distri- 
bution in the tail and trailing side of the disk. We find that 
the timescale for the HI tail to move from the location of the 
current UV tail to the current HI position is between 10-300 
Myr, depending on whether dense clouds decouple from the 
rest of the ISM. If the entire ISM stays together, then the star 
formation observed in the tail HII regions must have been ini- 
tiated when HI was located where the HII region is now. In 
this case the timescale for the HI tail to move is similar to an 
HII region lifetime, or <10 Myr. On the other hand, if some 
dense star-forming clouds decouple from the rest of the ISM 
(as seen in some other galaxies, e.g. Crowl et al. 2005), but 
most of the UV-bright stars form within the HI component of 
the ISM, then the star formation observed in the tail HII re- 
gions could have been initiated long after the lower density HI 
moved downstream. In this case the timescale for the HI tail to 
move is similar to the time since star formation stopped in the 
UV tail, or up to 200-300 Myr ago, based on optical-UV col- 
ors and stellar population models. Finally, if all star formation 
occurs in (long-lived) dense star-forming clouds which decou- 
pled from the HI (which we consider unlikely), the timescale 
could be even longer than 300 Myr. 

Other effects may be relevant for producing a UV-HI off- 
set, although further models are necessary to gauge the extent 
to which the tail angle might be affected. Turbulent motions 
at the ICM-ISM interface can cause the tail to shift positions, 
which could contribute to the HI tail's displacement from the 
UV/Ha tail. Judging from the hydrodynamic simulations of 
Roediger & Bruggen (2006, 2007), turbulent motions as the 
ICM strips away the ISM gas can easily cause a — 15° change 
in the angle between the denser parts of the tail and the ma- 
jor axis over timescales as short as —10 Myr. These short 
timescales are just within the lower limit of the amount of 
time NGC 4330's tail may have taken to shift positions. Al- 
ternatively, a combination of a non-linear star formation law 
plus projection effects can produce offsets between the peaks 
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of gas and star formation as observed in 2D maps (Vollmer et 
al. 2008), but the large-scale, systematic offset we observe in 
NGC 4330 makes this an unlikely explanation. Finally, rapid 
changes in ram pressure in a dynamic, structured ICM may af- 
fect the morphology of galaxies undergoing stripping. There 
is some evidence that shocks or bulk motions in the Virgo 
Cluster ICM can increase the ram pressure strength beyond 
what is expected from a smooth, static model (observations of 
Kenney et al. 2004; simulations of Tonnesen & Bryan 2008). 
More detailed modelling will be required to test whether the 
HI-UV offset implies slow or rapid changes in ram pressure. 

5.5. Extmplanar Star Formation 

Under what circumstances does star formation take place in 
the material stripped from galaxies, and where does it happen? 
Star formation is not sufficiently well-understood to be able to 
rely on models to predict where in the stripped gas it will take 
place. In NGC 4330, for at least some of the extraplanar UV 
regions, star formation must have occurred after their parent 
gas clouds left the disk. For instance. Region 5 (see Figurel23] 
for UV colors and Figure|2T|for location) has a stellar popula- 
tion age of ^10 Myr, and is located ^9 kpc from the galaxy's 
major axis. In order for star formation to have started while 
the cloud was still in the disk, the star-forming cloud would 
have had to move 9 kpc in 10 Myr, which translates to an 
average velocity of ^900 km s^' . Region 3 has a similar stel- 
lar population age and would require a velocity of ^600 km 
s^'. While these velocities are comparable to the ICM wind 
speed, they are probably much larger than the true velocities 
of the gas clouds, since clouds are not immediately acceler- 
ated to the ICM wind speed, and not all gas pushed out of the 
disk will achieve escape velocity. The fact that nearly all of 
the extraplanar star formation regions are associated with HI, 
which has velocities within 100 km s^' of the disk velocities 
(see Section l4~4l and FigurefTTT). is evidence that these regions 
are bound to the disk and have much lower velocities. It is 
therefore much more likely that at least some of the extrapla- 
nar star formation is initiated in extraplanar regions, and not 
in the disk. 

5.6. New Stellar Components Formed in Ram Pressure 
Stripped Gas 

In NGC 4330 we observe recent star formation that has oc- 
cured outside the thin disk due to ram pressure. Assuming 
that these stars formed from gas with roughly the same kine- 
matics as the extraplanar HI, most of these newly formed stars 
are close enough to the galaxy that they are likely bound to it, 
and will form new stellar components within the galaxy. On 
the trailing side, many UV-bright stars in the tail are displaced 
from the major axis by as much as 3.5 kpc. Over a few dynam- 
ical times, the stars that are currently UV-bright will spread 
out azimuthally and vertically to form a thick disk, with a ra- 
dial extent of > 12 kpc (the radius of the stars currently at 
the tip of the tail). Stars that form further from the disk, out 
along the tail or in other stripped gas, and remain bound to 
the galaxy, will create new halo components. In some cases 
the distinction between thick disk and halo will be arbitrary. 
The Ha upturn feature will also create a thick disk compo- 
nent, although possibly a much narrower one in radial extent 
than the one formed by the broad base of the tail on the trail- 
ing side. An upturn feature may propagate radially inward as 
stripping advances, injecting young stars into a thick disk or 
halo at progressively smaller radii. 



In contrast to thick disks made by mergers, which are older 
than their embedded thin disks the thick disks made by ram 
pressure stripping would be younger than their embedded thin 
disks. They are also expected to have a narrow range of stel- 
lar ages at any location. The vertical height and age of the 
stars formed through ram pressure are expected to vary with 
radius. Since galaxies are stripped from the outside in, in gen- 
eral younger thick disk or halo stars are expected closer to the 
galaxy center. The example of NGC 4330 demonstrates that 
ram pressure stripping can produce stellar thick disk and halo 
components with distinctive morphologies and age distribu- 
tions. Identifying such stellar components hold the potential 
for revealing the history of past ram pressure interactions in 
galaxies. 

6. CONCLUSIONS 

NGC 4330 provides an outstanding opportunity to study the 
effects of active ram pressure stripping on a nearby galaxy. 
The wind angle, evolutionary stage, and galaxy orientation 
make it possible to clearly observe distinct behavior at the 
leading edge (the upturn) and the trailing side (the tail) of 
the ICM-ISM interaction. Its stellar disk is undisturbed, with 
no warps, bridges, strong asymmetries, or other evidence of 
gravitational interaction. Ram pressure stripping is the only 
plausible explanation for its radially truncated, one-sided ex- 
traplanar HI distribution and undisturbed old stellar disk. 

We conclude by describing the morphological features re- 
lated to stripping, how we think they were created, and the de- 
rived evolutionary timescales for changes in the galaxy over 
the past ^400 Myr. Finally we provide a self-consistent 
description of the timeline of the galaxy's stripping history 
which takes into account these morphological features and 
their evolutionary timescales. 

1. Morphological Features Related to Stripping and 
Derived Stripping Interaction Parameters: 

• ISM truncation along the major axis: The galaxy's 
ISM has been stripped to well within the undis- 
turbed old stellar disk. Its dust extinction and 
HI are truncated at roughly the same radii, ^ 
50% — 65% of R25 on both the leading and trailing 
sides of the galaxy. 

• Extraplanar HI: The galaxy's HI distribution is 
strongly asymmetric as a result of having been re- 
moved from the disk by the ICM - 32% of the 
emission is extraplanar on the side we identify as 
downwind of the major axis. This includes 10% 
of the galaxy's HI in a long tail. 

• UV major-axis asymmetry: Recent star forma- 
tion as traced by UV is quite asymmetric about 
the galaxy's nucleus - the UV emission strength 
drops rapidly at 0.5 R25 on the upwind side, but 
gradually tapers off until R25 on the downwind 
side. The UV emission extends radially further 
on the upwind side than other tracers of the ISM 
and star formation - ongoing star formation as 
traced by Ha is truncated at ^ Q.5R25, roughly 
the same radii as the other ISM tracers, and is 
relatively symmetric about the galaxy's nucleus. 
The different extents of Ha and UV emission in- 
dicate that star formation extended to larger radii 
when the UV-bright stellar populations formed a 
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few hundred Myr ago, but was much more radi- 
ally truncated when the Ha-bright regions formed 
<10 Myr ago. Based on the UV color gradient in 
the outer disk on the galaxy's leading side, we es- 
timate that it has taken 200-400 Myr to strip the 
disk from >8 kpc to 5 kpc. However, there is 
no color/age gradient on the trailing side, so the 
effects of stripping have probably been more uni- 
form and recent on the trailing side. 

• The "upturn": At the leading edge of the ICM- 
ISM interaction, the galaxy's Ha emission veers 
abruptly out of the disk midplane, and a large, 
significantly obscuring dust cloud extends beyond 
the Ha upturn emission at a projected angle al- 
most perpendicular to the disk. This striking "up- 
turn" feature is near the HI truncation boundary, 
and the UV emission extends radially beyond it, 
suggesting that the ISM in the disk beyond the 
upturn was removed only recently. The upturn's 
Ha morphology indicates that most of the ISM in 
this region has been pushed out of the disk by ram 
pressure. The galaxy also displays a strong, local- 
ized radio deficit region near the upturn, charac- 
teristic of galaxies experiencing strong ram pres- 
sure. 

• The tail: NGC 4330 has a long tail on its trail- 
ing side which is visible in HI, UV, and Ha. The 
HI tail is offset downwind from the the UV/Ha 
tail, suggesting that much of the ISM has moved 
from the location of the UV tail to the current lo- 
cation of the HI tail. There is significant UV emis- 
sion from the tail but the Ha is sparse, indicating 
that there was significant star formation in the past 
'^10 - 200 Myr but there has been little in the past 
^5 Myr. We estimate that it has taken the ISM tail 
- 1 - 300 Myr to move downwind from the posi- xhis work was partially supported by the National Science 
tion of the UV-bright tail to the current position of Foundation under grant number 0607643 to Columbia Uni- 
the HI tail. This compai-es with the galaxy's rota- ygj-sity. GALEX (Galaxy Evolution Explorer) is a NASA 
tional period of -400 Myr at a radius of 5 kpc. Small Explorer, launched in April 2003. We gratefully ac- 

• Extraplanar UV-bright regions: We detect 9 re- knowledge NASAOs support for construction, operation, and 
gions downwind of the galaxy but outside the tail science analysis for the GALEX mission, developed in coop- 
that are bright in the U V, but for the most part lack eration with the Centre National d' Etudes Spatiales of France 
Ha emission. A comparison of their UV/optical and the Korean Ministry of Science and Technology. 

APPENDIX 

APPENDIX: COMPLICATING EFFECTS IN UV - OPTICAL COLOR INTERPRETATION 

In this section we discuss corrections for the three effects identified in Section 14.81 which complicate the interpretation of 
the UV-optical colors used to derive stellar population ages. Optical dust extinction is observed on the trailing side, which 
likely reddens the observed colors and cause the stellar populations to appear older than they actually are. The composite 
stellar population along the line of sight is a consequence of observing a nearly edge-on galaxy - within a single aperture, one 
measures stars at many galactocentric radii. If there is a stellar population gradient (as there is on the leading side of NGC 4330), 
one observes a composite stellar population, and it may be difficult to uniquely determine the star formation history as a function 
of galactocentric radius. Contamination from stars in the outer disk leads to redder or "older" colors being measured at a given 
radius. Finally, asymmetrical quenching times arise because stars now at the leading side were not always at the leading side, 
and may not have formed on the leading side, because the galaxy rotates. For example, any stars 150 - 300 Myr old which are 
now on the leading side between 5-8 kpc formed on the trailing side, since 200 Myr is half of the galaxy's rotational period. 
Some of the stars now on the leading side reflect the conditions on the trailing side in the past. 

1 . Dust Extinction. 

Figure [TSl includes a reddening vector calculated for AFuv=i, which corresponds to E(B-V)=0.035 using the wavelength- 
dependent dust extinction formulae of Calzetti (2001). This is a modest amount of reddening - for a location near R25 of 
a non-stripped late-type spiral, Muiioz-Mateos et al.(2009) find that Afuv=2- In particular, there is neither optical dust 



colors with stellar population models indicates 
that most of the regions have ages < 100 Myr, 
though the oldest is —350 Myr. Based on the ages 
and distances of the most distant regions from the 
galactic plane and the kinematics of the associated 
HI, it is likely that star formation in these regions 
was initiated after the gas left the disk. 

• The projected ICM wind angle: We discuss three 
methods to constrain the projected ICM wind an- 
gle from observations: the morphology of the 
large scale stripped gas tail, the orientation of 
elongated small-scale dust features, and the shape 
of the radio deficit region. We give quantitative 
estimates for all three techniques in NGC 4330, 
and constrain its disk-ICM wind angle to be 30- 
70°. 

2. Stripping History: The observations are consistent 
with the following scenario. Ram pressure increased 
rapidly — 200 — 400 Myr ago, stripping gas from the 
disk on the leading edge side from >8 to 5 kpc. The 
disk was stripped asymmetrically, and for a while the 
gas disk was stripped more deeply on the leading side 
than the trailing side. After — 1 rotation period, (—400 
Myr) at this higher pressure, the gas disk has become 
more symmetric, as currently observed. Within the last 
—400 Myr, as the gas disk has gone from highly asym- 
metric to fairly symmetric, gas on the trailing side has 
left the disk to form the base of the tail, and the tail 
has continued to move from the disk. It is not yet clear 
whether smooth or abrupt changes in ram pressure are 
responsible for this galaxy's recent, dramatic evolution. 
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Face-on exponential disk 




Fig. 24. — This figure illustates that when we observe an edge-on galaxy, we observe not only light from the annulus (green circles) at the proj ected galacto- 
centric radius (blue), but also hght from larger galactocentric radii that happens to be along the line of sight. For the apertures measured in Section l4!8l we must 
model the light distribution in the outer disk in order to determine how much of the observed flux is actually coming from the radius we are interested in. 

extinction nor mid-far IR dust emission on the leading side beyond 5 kpc (Abramson et al. 2009). NGC 4330 has had its 
outer disk stripped and now lacks any ISM tracers on the leading side. As a result of stripping, the outer disk probably has 
much less dust extinction than a typical spiral. 

2. Composite Stellar Population Along the Line of Sight. 

Because we view the galaxy nearly edge-on, each aperture measured in Section 14.81 contains not only light from the 
projected galactocentric radius, but also light from further out in the disk which happens to be along our line of sight 
(Figure l24li. To correct for this effect, we model the radial light distribution of the outer disk and subtract the amount of 
light in each band which comes from a radial bin further out in the disk than the one we are interested. We used Galfit 
(Peng et al. 2002) to fit an exponential disk to the galaxy's light distribution at each wavelegth, with separate fits to the 
upturn and tail sides. By doing this, we can estimate what fraction of the flux in each aperture is actually coming from the 
radial bin that we are interested in. The different scale lengths of the UV and r-band disks result in different correction 
factors for each wavelength, changing the FUV-NUV and NUV-r colors. 

We also tried fitting a power-law function, and a combination of exponential and power-law functions, and found that every 
function that provided a very good match to the data produced similar correction factors. However, if the actual disk is 
non-axisymmetric, which is likely, the correction factors could be somewhat different. There are several possible ways to 
fit each light distribution, and the derived functions are most affected by which portion of the disk's radial light distribution 
is fit. Although it is difficult to uniquely determine the star formation history at a given galactocentric radius, we can use 
this technique to estimate how much the observed colors might change if we were able to remove the light from the outer 
disk, which allows us to estimate the range of possible quenching times for the outer disk. 

The range of possible values resulting from the outer disk correction are given by the colored boxes in Figure [18] Without 
applying the correction, the derived stellar population quenching times (Figure[^ filled symbols) are quite old (250 - 400 
Myr), even those at 5 kpc, just beyond the upturn structure at the leading edge. Three Ha-bright regions at the edge of the 
upturn structure have uncorrected quenching times of between 50 and 200 Myr, although the lifetime of an HII region is 
generally closer to ~3 Myr. After applying the correction (the parameter space within the open boxes in Figure [IST i, the 
point closest to the upturn at 5 kpc may have a quenching time of <C50 Myr, which is much more realistic. The correction 
introduces additional error which is hard to estimate, but it provides a general idea of how much the outer disk may redden 
the colors in the apertures we observe. 

3. Asymmetrical Quenching Times. 

The quenching time as a function of position does not have a simple relationship with the star formation history as a 
function of position because stars move. Stars now at the leading side were not always at the leading side, and may not 
have formed on the leading side, because the galaxy rotates. For example, some of the stars now on the leading side 
reflect the conditions on the trailing side in the past: any stars 150 - 300 Myr old which are now on the leading side at 
radii between 5-8 kpc formed on the trailing side, since 200 Myr is half of the galaxy's rotational period. This leads to 
additional uncertainty in the stripping rates we derive. NGC 4330's maximum HI rotational speed is 130 km/s, so for a 
radius of 5 - 8 kpc, the rotational period is 230 - 380 Myr. This is comparable to the quenching time gradient on what 
currently constitutes the leading side, ^300 Myr, so the stellar populations we observe at the leading side of the galaxy 
were not simply quenched in situ at their current location. 

In order to approximate a gradual quenching of star formation, as might happen if a given stellar population suffered 
somewhat suppressed star formation when the ISM was partially removed at the start of stripping, then kept evolving and 
was fully quenched on the leading side, we have used models with two episodes of quenching. First, the star formation 
rate was decreased by 50% or 90% either 200 or 500 Myr ago. Then the star formation was completely quenched at a 
range of more recent times - the results are the colored lines in Figures |25] The highest log(FUV/NUV) value for each line 
represents a quench of the remaining star formation 50 Myr ago, with the second quench happening progressively longer 
ago. For example, the solid orange line represents a 90% quench 200 Myr ago, and if the remaining 10% of star formation 
was quenched 50 Myr ago, the colors are very similar to a full quench 200 Myr ago. 
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Fig. 25. — Same as Figure [Tsl but with colored lines indicating partially-quenched star formation histories. The solid lines indicate stellar populations whose 
star formation rates were quenched by 90% either 500 Myr (blue) or 200 Myr (orange) ago, and had the remaining star formation quenched more recently. The 
dashed lines indicate a 50% quench at 500 Myr (blue) and 200 Myr (orange) ago. 

Interestingly, the two-quench models are not a good fit to the observed UV/optical colors - they are a worse fit than the 
quench plus burst model. 
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